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SLQIDES,  TWO-DIMENSIONAL  B(»IES  HAVINQ  . 
MINIMUM  TOTAL  DRAG  AT  HYPERSONIC  SPEEDS 

by 

ANGELO  MIEI£^*^  and  ROBERT  E.  HIITCHARD^**^ 


SUMMARY 

This  paper  considers  the  problem  of  minimizing  the  total  drag  (sum 
of  the  pressure  drag  and  the  friction  drag)  of  a  slender,  two-dimensi.onal, 
symmetric  body  at  zero  angle  of  attack  in  hypersonic  flow  under  the  as¬ 
sumption  that  the  distribution  of  pressure  coefficients  is  Newtoniaui  and 
that  the  friction  coefficient  is  constauit*  After  the  condition  that  the 
pressure  coefficient  must  be  nonnegative  is  accoimted  for,  the  minimal 
problem  is  solved  for  arbitrary  conditions  Imposed  on  the  thickness,  the 
enclosed  eorea,  and  the  moment  of  inertia  of  the  contour  under  the  assump¬ 
tion  that  the  length  is  free.  It  is  shown  that,  if  convenient  dimension¬ 
less  coordinates  are  employed  (that  is,  if  the  abscissa  and  the  ordinate 
are  normalized  with  respect  to  the  length  and  the  semi  thickness),  the  to-  . 
tality  of  extremal  arcs  is  composed  of  a  two-parameter  family  of  solutions. 
It  is  also  shown  that  each  extremal  arc  involves  at  most  one  comer  point 
and,  hence,  two  subarcs:  one  of  these  is  characterized  by  a  positive  pres¬ 
sure  coefficient  and  is  called  the  regular  shape;  the  other  is  characterized 
by  a  zero  pressure  coefficient  and  is  called  the  zero-slope  shape,  Thus, 
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two  classos  of  bodies  can  be  identified:  (I)  bodies  composed  of  a  regular 
shape  onljt  and  (II)  bodies  composed  of  a  regular  shape  followed  by  a  con¬ 
stant  thickness  contour* 

Particular  attention  is  devoted  to  solutions  for  which  either  one  or 
two  of  the  quantities  under  consideration  are  prescribed*  If  only  one 
quantity  Is  given  (the  thicknessf  the  enclosed  area,  or  the  moment  of  inertia 
of  the  contour)  the  extremal  arc  is  a  single  curve  of  class  I  regardless  of 
the  friction  coefficient*  On  the  other  hand,  if  two  geometric  quantities 
are  given  (the  thickness  and  the  enclosed  eurea,  the  thickness  smd  the  moment 
of  inertia  of  the  contour,  or  the  enclosed  area  and  the  moment  of  inertia 
of  the  contour),  a  one-parameter  family  of  extremal  arcs  exists;  the  pa¬ 
rameter,  called  the  friction  parameter,  is  proportional  to  tL.  cubic  root 
of  the  friction  coefficient  and  is  related  to  the  quantities  which  are  pre¬ 
scribed*  Depending  on  the  value  of  the  friction  parameter,  two  distinct  be¬ 
haviors  are  possible*  If  the  friction  parameter  is  subcrltical  (smaller  than 
a  certain  critical  value),  the  solution  is  of  class  I;  if  the  friction  pa¬ 
rameter  is  supercritical  (larger  than  a  certain  critical  value),  the  solution 
is  of  class  II  with  the  transition  point  from  the  regular  shape  to  the  con¬ 
stant  thickness  contour  shifting  forward  as  the  friction  parameter  increases* 
For  all  of  the  cases  considered,  analytical  expressions  are  derived  for  the 
optimum  shapes,  the  thickness  ratios,  and  the  drag  coefficients* 
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1#  INraODUCTIQN 

The  problem  of  ninimizing  the  drag  of  slender,  two-dimensional  bodies 
in  hypersonic  flow  has  attracted  considerable  attention  in  recent  times* 

With  particular  regard  to  the  pressure  drag,  general  solutions  have  been 
obtained  in  Ref*  1  under  the  assumption  that  the  pressure  distribution  is 
Newtonian  and  that,  among  the  geometrical  quantities  being  considered  (the 
thickness,  the  length,  the  enclosed  area,  and  the  moment  of  inertia  of  the 
contour),  two  are  prescribed  and  the  remaining  two  are  free*  These  solu¬ 
tions  have  been  extended  in  Ref*  2  to  cover  the  case  where  three  of  these 
quantities  are  given  and  only  one  is  free* 

While  the  investigations  of  Refs*  1  and  2  neglected  the  friction'  drag, 
it  should  be  noted  that  there  exist  practical  values  of  the  thickness  ratio 
for  which  the  friction  drag  may  have  the  same  order  of  magnitude  as  the  pres¬ 
sure  drag*  Therefore,  it  is  of  Interest  to  reinvestigate  the  problem  of  the 
optimum  slender  shape  from  the  point  of  view  of  minimizing  the  total  fore- 
body  drag  (sum  of  the  pressure  drag  and  the  friction  drag)  for  any  number 
of  conditions  imposed  on  the  thickness  d,  the  length  /,  the  enclosed  area  A, 
and  the  moment  of  inertia  of  the  contour  M*  This  is  the  problem  considered 
in  the  present  report  in  connection  with  the  following  assumptions: 

(a)  the  body  is  slender;  (b)  the  distribution  of  presswe  coefficients  is 
Newtonian;  auid  (c)  the  friction  coefficient  is  constant  along  the  contour* 

The  corresponding  axisyonetric  problem  is  analyzed  in  Ref*  3  for  any  number 
of  conditions  Imposed  on  the  diameter,  the  length,  the  wetted  area,  and  the 
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2*  MINIMIIM  DRAG  PROBlSi 

CoDBlder  a  two-dlmenBloiialy  BTametric  airfoil  at  zwo  angle  of  attack 

in  a  hyperaonlc  flow,  and  denote  by  x  a  coordinate  in  the  flow  direction, 

y  a  normal  coordinate,  and  y  the  derivative  dy/dx*  Uhder  the  alender  body 
2 

approximation  y  «  1,  the  assumed  Newtonian  distribution  of  preseure  coef- 

2 

flcients  simplifies  to  C  «  2y  «  Consequently,  the  drag  per  unit  span  of 

P 

that  portion  of  the  body  which  is  included  between  stations  0  and  x  is  given 

by 

DW  .  dx  (1) 

where  is  the  friction  coefficient,  assumed  constanto  The  corresponding 
values  for  the  area  enclosed  by  the  contour  and  the  moment  of  inertia  of 
the  contour  are  given  by 

A(x)  ■  2  r  ydx  ,  M(x)  »  2  T  y^dx  (2) 

Jq  Jo 


After  the  definitions 


t.Hf,  y.SSf.  (3) 

aure  introduced,  differentiation  of  both  sides  of  £k|S»  (1)  and  (2)  with  re«» 
spect  to  the  independent  variable  leads  to  the  following  differential  con¬ 


straints  s 
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P  -  y  =  0  (4) 

.  2  ^ 
y  •  J  a  0 


Since  the  requirement  that  the  eJLope  be  nonnegative  everywhere  can  be  ex- 
preesed  as 


y  - 


0 


(5) 


where  p  denotes  a  real  variable ^  the  differential  eystem  composed  of  Eqs*  (4) 
and  (5)  involves  one  independent  variable  (x),  five  dependent  variables 
(yp  Oft  09  Y,  p)9  and  one  degree  of  freedom*  In  this  connection 9  after  as** 
Burning  that 


«  er^  «  $i  «  Yi  ■  0  (6) 

and  that  some  9  but  not  all9  of  the  remaining  state  variables  are  given  at 
the  final  point9  one  can  formulate  the  minitniifn  drag  problem  as  follows: 

In  the  class  of  functions  y(x),  cKx),  B(x)*  v(x)*  p(x)  which  are  consistent 
with  the  differential  constraints  (4)  and  (5)  and  the  initial  conditions  (6)^ 
find  that  special  set  which  minimizes  the  difference  ^  »  Q^*  wnere 


Q  »  Of* 
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3*  WECESSAHT  COWDITIOKS 

Th«  previous  problem  is  of  the  Mayer  type  with  separated  end  oonditions* 
Consequentlyt  after  the  Lagrange  multipliers  throu^  Intredueed 

and  the  fundamental  function  Is  written  as  (Hefs.  3  and  6) 

y  ■  Xi  ^or  -  y^  -  +  XgCP  -  y)  +  Xj(y  -  y^)  +  \^r  “  P^)  (7) 

the  extremal  aro  is  described  by  the  following  Guler»Lagrange  equations! 


^(X4  -  -  -  Xg  -  eXjy 

1^-0 

Xg  -  0  (8) 

Xj  -  0 

0  -  X4P 

the  second,  third,  and  fourth  of  idilch  can  be  Integrated  to  give 

Xi  -  Cl  ,  Xg  -  Cg  ,  Xj  -  Cj  (9) 


where  through  are  constants*  Furthermore,  after  it  is  observed  that 
the  fundamental  function  does  not  contain  the  Independent  variable  expli¬ 
citly  and  after  Bq*  (8-3)  is  accounted  for,  the  following  first  integral  can 
be  established! 
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(10) 


where  C  Is  a  conetanto 

Comer  conditions#  As  the  fifth  Euler  equation  indicates,  the  extre¬ 
mal  arc  is  composed  of  the  subaros 

e  0  and/or  p  ■  0  (11) 

Along  the  former  subarcs,  called  regular  shapes »  the  pressure  coefficient 
is  always  positive  as  long  as  p  is  real*  Along  the  latter  subarcs,  called 
nero-slope  shapes,  the  pressure  coefficient  is  always  zero*  The  Junctions 
between  the  subau'cs  must  be  studied  with  the  aid  of  the  Erdmann-Weierstrass 
comer  conditions*  They  require  that  the  integration  constants  C^,  C^t  0 
have  the  same  value  for  each  of  the  subarcs  composing  the  extremal  arc  and 

&(P)  -  A(X4  -  -  0  (12) 


where  A(**«)  denotes  the  difference  between  quantities  evaluated  after  and 
before  the  corner  point*  A  mathematical  consequence  of  these  equations 
are  the  relationships 

AX4  «  Ay  -  0  (13) 

^*^Eq«  (12-1)  is  a  consequence  of  the  first  integral  (10)  and  the  con¬ 
tinuity  of  the  integration  constants*  , 
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whlohf  if  ooablaed  with  Eqa*  (11),  iflply  that 

Xjj  ■  y  ■  0  (14) 

before  and  after  a  corner  point* 

End  conditions*  The  end  conditions  are  partly  of  the  fixed  end-point 
type  and  partly  of  the  natural  type*  %e  latter  nnist  be  deterained  froa 
the  transversality  condition 

j^-  Cdx  +  (Cj^  4-  1)  d«  C2dP  ♦  Cjdy  *  i\  ~  dyj^  «  0  (15) 

which  must  be  satisfied  for  every  system  of  differentials  consistent  with 
the  prescribed  end  conditions;  in  particular,  it  implies  that  1* 

If  the  length  is  free,  the  transversal! ty  condition  yields  C  «  0*  On 
the  other  hand,  if  either  the  enclosed  area  or  the  moment  of  inertia  of  the 
contour  are  free,  the  transversal! ty  condition  leads  to  ■  0  and  «  0, 
respectively*  Finally,  if  the  diameter  is  free,  the  transversality  condition 
leads  to 


(X4  ♦  3y^),  -  0 


(16) 


which,  if  combined  with  the  Euler-Iagrange  equation  (8-^),  iiq>lies  that 


\t  -  ^f  •  0 


(17) 
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Conaaquentlx,  if  d  »  denotes  the  dleneteri  the  first  integral  (10) 
yields  the  additional  relationship 

A 

C+y-C2^-Cj-jj-»0  (18) 

At  this  point,  it  is  convenient  to  separate  the  discussion  into  two 
basic  problems:  problems  where  the  length  is  given  and  problems  where  the 
length  is  free*  As  Eq*  (1)  shows,  problems  of  the  first  kind  are  char¬ 
acterized  by  the  fact  that  the  friction  drag  is  independent  of  the  shape  so 
that  the  contour  which  minimizes  the  total  drag  is  identical  with  that  which 
minimizes  the  pressure  drag*  Since  shapes  of  minimum  pressure  drag  have 
been  fully  discussed  in  Refs*  1  and  2,  these  problems  are  not  considered 
here*  Thus,  the  analysis  is  restricted  to  problems  of  the  second  kind,  in 
which  the  length  is  free*  The  class  of  problems  in  which  the  length  is  free 
contains  several  subclasses  which  depend  on  the  number  of  quantities  that  are 
specified*  Among  these  subclasses,  the  following  are  considered  here: 
problems  in  which  one  geometric  quantity  is  given  and  problems  in  which  two 
geometric  quantities  are  given*  For  these  problems,  simple  manipulations 
lead  to  the  results  which  are  susmarized  in  Table  1  where  two  types  of  re¬ 
latione  are  indicated:  those  obtained  from  the  transversality  condition  and 
those  obtained  by  combining  the  results  of  the  transversality  condition  with 
the  Euler-Iagrange  equation  (8-^)  and  the  first  integral  (10)* 

Legendre-Clebsch  condition*  !nie  Legendre-Clebech  condition  indicates 
that  the  drag  is  a  minimum  if  the  following  inequalities  are  satisfied  every¬ 
where  along  the  extremal  arc: 
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j  >  0  ,  along  tha  regular  shapa 

<  0  •  along  tha  naro-alopa  ahapa 


(19) 


Switching  function#  From  tha  previoua  diacuaaiont  it  appear  a  that  tha 
lagranga  multipliar  X|^  plaja  an  important  role  in  datarmining  tha  compoaition 
of  tha  axtranal  arc*  If  tha  tarminologj  of  control  theory  ia  aaplojadt  thia 
multipliar  can  be  called  tha  switching  function!  its  properties  are  as  fol» 
lowas 

X|^  «  0  9  along  tha  regular  shape 

X|^  <0  9  along  the  nero-alopa  ahapa  (20) 

I4  *  0  9  at  a  cornar  point 
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4.  OEOHETRT  OF  THE  EXTREMAL  ABC 

In  the  previous  sectiont  the  necessary  conditions  to  be  satisfied  by 
the  extremal  arc  have  been  stated*  In  this  sectlont  several  general  conse¬ 
quences  of  these  equations  are  derived  referring,  for  the  sake  of  brevity, 
to  the  fflinimum  drag  problem  (C^  «  -  1)  with  the  length  unspecified  (C  •  0). 
In  order  to  facilitate  the  analysis,  the  following  dimensionless  coordinates 
are  introduced: 


?  -  x/i  ,  T|  -  2y/d 

together  with  the  definitions 

Kg  » 


(21) 


(22) 


With  these  coordinates,  the  first  integral  (10)  reduces  to  the  form 
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1 


(23) 


where  t  »  d/i  denotes  the  thickness  ratio  and  1)  the  derivative  dT/d^« 

Basic  inequatlity*  The  application  of  the  above  first  integral  at 
the  end  points  of  the  extremal  arc  indicates  that  the  terminal  values  of  the 
slope  are  given  by 
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(24) 


From  the  first  equatioiii  it  Is  apparent  that  the  optiaum  body  is  sharp- 
nosed*  Since  the  final  slope  (and* hence,  the  final  pressure  coefficient) 
oust  be  nonnegative,  one  deduces  ftom  the  second  equation  that  the  follow¬ 
ing  basic  inequality  must  be  satisfiedi 

1  -  K2  -  K j  >  0  (25) 

Switching  function*  Since  each  extremal  arc  may  involve  more  than  one 
subaro,  it  is  of  paramount  importance  to  calculate  the  distribution  of  the 
switching  function;  in  nondimensional  form,  this  function  can  be  defined  as 

a  -  ^  (26) 

For  the  regular  shape,  it  is  known  that  a  >  0*  For  the  nero-slope  shape,  it 
is  known  that  f|  ■  0  and  1]  «  const*  Consequently,  the  Euler-Iagrange  equation 
(8-1)  reduces  to 


a  -  -  K2  -  2KjT) 


(27) 


which,  in  the  light  of  the  initial  conditions  (20-5),  admits  the  particular 
integral 
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a  -  -  (K2  +  2KjT^)  (?  -  ?^)  (28) 

1a  which  the  subscript  c  refers  to  a  corner  point# 

Sequence  of  subarcs#  The  next  step  is  to  determine  the  appropriate 
sequence  of  subarcs#  First,  it  is  observed  that  the  extremal  arc  cannot 
start  with  the  zero-slope  shape  T]  ■  0:  the  equation  of  this  shape  would 
be  incompatible  with  the  first  integral  (23) •  Thus,  if  it  is  assumed  that 
the  extremal  arc  starts  with  a  regular  shape,  the  next  questions  to  be  in¬ 
vestigated  are:  (a)  whether  a  transition  to  a  zero-6loi>e  shape  is  possible; 
and  (b)  if  so,  whether  a  further  transition  back  to  a  regular  shape  may  occur* 
Concerning  the  first  question,  the  comer  condition  (l4-2)  and  the  first  in¬ 
tegral  (23)  show  that  the  transition  from  the  regular  shape  to  the  zero-slope 
shape  is  possible  if  the  following  relationship  is  satisfied: 

1  -  K2T^  -  KjT]^  =  0  (29) 

With  regard  to  the  second  question,  Eq*  (28)  shows  that,  since  the  switching 
function  varies  linearly  with  the  abscissa  along  the  zero-slope  shape,  it 
can  only  vanish  at  one  point:  the  corner  point  between  the  regular  shape 
and  the  zero-slope  shape©  Thus,  no  regular  shape  may  follow  the  zero-slope 
shape  so  that  the  equation  of  the  latter  is  T)  s  1«  Furthermore,  because 
of  £q©  (29)  and  the  properties  of  the  switchdLng  function,  the  presence  of 
a  zero-slope  shape  requires  that 

1  -  K2  -  Kj  ■  0 


♦  2Kj  ^  0 


(30) 


I4 


5ino#  no  noro  tham  one  comer  point  and  two  subaroa  may  exiatf  the  totality 
of  extremal  area  conalsts  of  two  claaaea  of  bodieai  (I)  bodlea  conpoaed  of 
a  regular  shape  only,  and  (II)  bodlea  compoaed  of  a  regular  ahape  followed 
by  a  constant  thlckneaa  contour*  Iheae  bodies  are  repreaented  symbolically 

by 


Class  It  0*0 

Class  II  t  o*0-*T)»l 


(31) 


Faunlly  of  solutions*  Since  the  most  general  ^pe  of  extremal  arc  is 
of  class  III  its  geometry  can  be  described  by  the  equations 

^  (1  -  K2T1  -  dll 

(1  -  K^T)  -  dll  (32) 

5^  <  ?  <  1  ,  T1  -  1 


0  <  ?  < 


I 


L 


where  I  denotes  the  abscissa  of  the  transition  point*  Bodies  of  class  I 
can  be  obtained  from  bodies  of  class  II  by  means  of  the  formal  substitution 
^  ■  1;  it  ahould  be  no  ted  t  however  •  that  the  comer  condition  need  not  be 
aatiafied  at  this  special  point*  In  a  functional  form^  Eqs*  (3^)  can  be  re-¬ 
written  as 


iK§.  5ct  h*  V 


(33) 
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60  that,  after  this  relationship  is  combined  with  either  of  the  relations 

Clasa  I  :  ^  «  1 

(54) 

Class  II  t  1  -  •  Kj  »  0 

it  is  seen  that  a  two-parameter  family  of  optimum  bodies  exists* 

For  particular  types  of  boundary  conditions,  considerable  simplifica¬ 
tions  are  possible*  Thus,  if  the  enclosed  area  is  free  (K^  *  0)  or  the 
moment  of  inertia  of  the  contour  is  free  (K^  *  0),  the  number  of  independ¬ 
ent  parameters  is  reduced  by  one*  An  analogous  remsork  holds  for  the  case 
where  the  thickness  is  free,  since  K2  Kj  *  1*  In  conclusion,  the  number 
of  independent  parameters  governing  the  solution  depends  on  the  number  of 
geometric  quantities  other  than  the  length  which  are  prescribed*  If  three 
quantities  are  prescribed,  the  problem  admits  a  two-parameter  family  of  so¬ 
lutions*  If  two  quantities  are  prescribed,  the  problem  admits  a  one- 
parameter  family  of  solutions*  Finally,  if  only  one  quantity  is  prescribed, 
the  problem  admits  a  zero-parameter  family  of  solutions,  that  is,  the  geometry 
of  the  extremal  arc  in  the  ^l^plane  consists  of  a  single  curve  regardless  of 
the  value  of  the  friction  coefficient*  In  connection  with  problems  where 
either  one  or  two  geometric  quantities  are  prescribed,  the  set  of  dimension¬ 
less  boundary  conditions  is  indicated  in  Table  2  along  with  the  dimensionless 
switching  function  at  the  final  point,  the  slope  of  the  extremal  arc  at  the 
final  point,  and  the  number  of  independent  parameters  governing  the  solution* 
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5.  SOLOTIOW  OF  IHE  BOUNDARY  VALUE  PROBLEM 

In  thin  sectloui  a  general  method  for  determining  the  unkaoima  ap¬ 
pearing  in  Eqs«  (32)  is  presented*  The  analysis  is  facilitated  if  sereral 
nondimensional  integrals  are  introduced*  If  the  cubic  root  of  both  sides 
of  £q*  (23)  is  extracted,  if  the  variables  are  separated,  and  if  an  inte¬ 
gration  over  the  regular  shape  is  performed,  the  following  result  is 
obtained! 


K2,  Kj) 


(55) 


tAtr* 


r  •  i/ 1 


(36) 


daaotas  tha  thicknasa  ratio  and  the  nondifflanaional  Intagral 


Id<5c»  *2.  K3) 


f  a  )"^  dTl 

^c  Jo  ^ 


(37) 


Furthermore,  by  simple  manipulations,  the  area  enclosed  by  the  airfoil 
and  the  moment  of  inertia  of  the  contour  become 


AT/d^  -  Kg,  Kj) 

2Hr/i^  “  ''2*  *3^ 


(38) 


%diere 
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V'o’  *2 


■  ■£ 


T>i?  -  1  -  5, 


(1  -  Tl)  (1  -  Kgl)  -  dT) 


c 

0 


y  (1  -  Kgn  -  Kjlf  )“^  dT) 


(39) 


I„(?c.  K^.  K,)  -j  D^d?  .  1  -  -0— - 

J  (1  -  K2T)  -  dT) 


Fbr  a  given  friction  coefficient |  the  eystein  composed  of  the  five  equations 
(3^)  through  (36)  and  (38)  involves  the  eight  quantities 

T*  i,  A,  Kj  (^) 

which  means  that  one  particular  optimuin  body  can  be  determined  if  three 
additional  relationships  are  specified*  For  the  boundary  conditions  con¬ 
sidered  in  !I^ble  2^  these  relationships  are  represented  by  amy  one  of  the 
following  sets: 


d  -  Const  a  0  , 
A  -  Const  a  0  I 
N  -  Const  a  0  9 
d  -  Const  a  0  • 
d  -  Const  a  0  , 
A  -  Const  a  0  , 


K2  a  1 
K2  a  0 

A  *  Const  a  0 
M  -  Const  a  0 
M  -  Const  a  0 


Kj  a  0 
K3.0 
K3.1 
Kj  »  0 


(41) 
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Drag  oo#fflolent»  Aft«r  tha  boundary  valua  problaa  haa  baan  aolTady 
tha  Dfxt  atap  ia  to  datarmlna  tha  drag  of  tha  optiaum  body*  Diia  drag  can 
ba  writtaa  aa 


whara  danotaa  tha  diaanaionlaaa  integral 
Dp 

hp  ■  <■>’ 

NoW|  if  tha  drag  coafficlent  la  rafarrad  to  the  frontal  area  at  x  ■  i  (that 
ia,  if  G|^  «  D/qd)^  the  following  relatlonehip  can  ba  readily  astabliohad 
batwaan  the  drag  coefficient |  tha  friction  coefficient  9  and  tha  thieknaea 
ratio I 


^  +  2  (44) 

Notice  thatf  if  both  aides  of  Eq*  (23)  are  multiplied  by  d^  and  integrated 
over  tha  entire  length  of  tha  extremal  arc,  tha  ralationahip 

4  ^ 

can  be  obtained*  Conaaquentlyi  after  Eqa*  (33)9  (^)9  and  (43)  are  eonbinedi 
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on*  deduct*  tbat 


T  -  T  -  ''2IA  -  W 


(46) 


Drag  ratlo<  Another  Intereetlng  quantity  characterizing  the  optljrai 
body  la  the  drag  ratio ,  that  Ibi  the  ratio  of  the  friction  drag  to  the  total 
drag*  Because  of  Iki*  (44 )«  this  quantity  ia  given  by 


4C, 


Sf  _  '"f 

4Cj  +  •?! 


(47) 


which,  in  th*  light  of  Eq.  (4$),  can  b*  rewritten  as 


’•'Df  2  - 

<=»  ■  3  -  V,  -  KjIm 


(48) 
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6.  PARTICULAR  CASES 

In  th«  previous  sections,  the  raininum  drag  problem  was  solved  in  general 
for  arbitrary  boundary  conditlonse  Here,  several  particular  eases  are  con¬ 
sidered,  and  the  associated  optimum  shapes  are  cadculated*  Two  classes  of 
problems  are  considered t  (a)  problems  in  which  only  one  geometric  quantity 
is  prescribed  and  (b)  problems  in  which  two  geometric  quantities  are  pre¬ 
scribed*  Problems  of  type  (a)  are  chax*acterized  by  solutions  of  class  I| 
that  is,  solutions  involving  a  regular  shape  only*  Problems  of  type  (b) 
are  characterized  by  solutions  of  either  class  I  (regular  shape)  or  class  II 
(regular  shape  followed  by  a  constant  thickness  contotir)  depending  on  whether 
the  friction  coefficient  is  smaller  or  larger  than  a  certain  critical  value* 

6*1*  Given  Thickness 

If  the  thickness  is  given  while  the  enclosed  area  and  the  moment  of 
inertia  of  the  contour  are  free,  the  transversality  condition  leads  to 

»  0*  Since  £q*  (30-1)  is  not  satisfied,  a  zero-slope  shape  cannot 
exist*  Hence,  the  extremal  arc  is  of  class  I,  that  is,  involves  a  regular 
shape  only* 

After  setting  ^  ■  1,  the  equation  of  the  regular  shape  (32-1)  can  be 
'C 

integrated  to  give  (Fig*  1) 


71-5  ('♦9) 

meaning  that  the  extremal  arc  is  a  wedge*  Since  I^  ■  1,  Eq*  (33)  yields 
the  following  value  for  the  optimum  thickness  ratio: 


21 


T  -  -  1.26  (50) 

Consequently,  Eq.  (46)  indicates  that  the  drag  coefficient  per  unit  thick* 
ness  ratio  squared  is  given  by 


T 


a  result  which,  because  of  Eq.  (48),  has  the  following  implicatioat  the 
friction  drag  of  the  optimum  body  is  two^thirds  of  the  total  drag# 

6.2.  Given  Enclosed  Area 

If  the  enclosed  area  is  given  while  the  thickness  and  the  moment  of 

inertia  of  the  contour  are  free,  the  transversality  condition  leads  to 

*■  If  *  0,  and  >  0.  Should  a  zero-slope  shape  exist,  the  s%dLtching 

function  would  be  zero  at  both  ends  of  this  subarc.  However,  because  of 

Eq.  (28),  this  is  only  possible  when  %  •  ^  •  Since  the  length  of  the  zero- 

X  c 

slope  shape  is  zero,  the  extremal  arc  is  of  class  I,  that  is,  includes  a 
regular  shape  only* 

For  ^  «  1,  the  equation  of  the  regular  shape  (32-1)  can  be  integrated 
c 

to  give  (Fig.  1) 


1)  -  1  -  (1  -  5)5/^  (52) 

Since  ■  3/2,  Eq*  (33)  yields  the  following  value  for  the  optimun  thickness 


ratios 
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T  ■  I  ^  -  0.84  ^  (53) 

Finally,  sine#  ■  3/5#  £q*  (^)  yields  the  following  ainiaua  drag  ooef- 
fioientt 

^  ^  ■  ^#05  (5^) 

T 

which,  in  the  light  of  Eq»  (48),  has  the  following  implication!  the  ftiotion 
drag  of  the  optimum  body  is  five^sixths  of  the  total  drag* 

6^3^  Given  Moment  of  Inertia  of  the  Contour 

If  the  moment  of  inertia  of  the  contour  is  given  while  the  thickness 
and  the  enclosed  area  are  free,  the  transversallty  condition  leads  to 

”  0,  «  1,  and  «  0*  Should  a  zero-slope  shape  exist,  the  switching 

function  would  be  zero  at  both  ends  of  this  subarc*  However,  because  of 
Eq*  (28),  this  Is  only  possible  when  Since  the  length  of  the  zero- 

slope  shape  is  zero,  the  extremal  arc  is  of  class  I,  that  is,  includes  a 
regular  shape  only* 

For  •  1,  one  can  integrate  the  equation  of  the  regular  shape  (32-1) 
to  obtain  (Fig«  1) 


? 


(55) 


where 
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f(T))  .  -  ~  ?•  r(«p,  k)  ♦  E(«p,  k)  - - ^  — -  (56) 

^  ^  1  -  gi  -  if 

and  where  F  and  £  denote  the  incomplete  elliptic  integrals  of  the  first 
and  second  kind,  respectively*  The  associated  argument  (p  and  parameter  k 
are  defined  as 


^1))  ■  arc  cos 


J3-1  +  Vl  - 

Tj-M-yrr?  ’ 


(57) 


Since  ■  3f(l),  Eq.  (35)  yields  the  following  value  for  the  optiaum  thick¬ 
ness  ratio I 


T  ■  jTfiy  “  0*97 


(58) 


Finally,  since  •  3/7i  the  linlmum  drag  coefficient  becomes 

^  f^(l)  -  2.78  (59) 

a  result  which,  because  of  £q«  (^),  has  the  following  implications  the 
friction  drag  of  the  optiaum  body  is  seven-^ninths  of  the  total  drag* 

6*4*  Given  Thickness  and  Enclosed  Area 

If  the  thickness  and  the  enclosed  area  are  prescribed  while  the  moment 
of  inertia  of  the  contour  is  free,  the  transversal! ty  condition  leads  to 
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«  0«  The  totality  of  extremal  arcs  is  represented  by  a  one-parameter 
family  of  solutions  of  either  class  I  (regular  shape)  or  class  II  (regular 
shape  followed  by  a  constant  thickness  contour)  depending  on  whether  the 
friction  coefficient  is  smaller  or  larger  than  a  certain  critical  value# 
The  representation  of  the  results  becomes  simple  and  immediate  if  a  thick¬ 
ness  parameter  and  a  friction  parameter  are  introduced*  These  parameters 
are  defined  by 

K  *  T  -4  t  Kf  ■  \ 

T  d  d 


and,  because  of  Bqs*  (35)  and  (38-1),  can  be  revrritten  as 


(61) 


Bodies  of  class  I.  These  bodies  consist  of  a  regular  shape  only  and 
are  obtained  for  f  «=  1  and  -  •  «  K-  s  1.  After  Eqs.  (32-1),  (46),  and  (6l) 

‘C  & 

are  employed,  the  optimum  shape,  the  thickness  parameter,  the  drag  coefficient, 
and  the  friction  paurameter  can  be  rewritten  as 


%  • 


1  -  (1  -  K^tl) 

1  -  (1  -  K2) 


2/3 


VT 


^  3  -  (3  +  2K^)  (1  - 

5*2  1  -  (1  - 


(62) 


(1  -  -  *2>  "  *2^^] 


K 


f 


-  (3  ♦  2K2)  (1  -  K2)^] 
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Elinlnatioa  of  the  paraneter  K2  from  these  equations  yields  the  funotlonal 
relationships 


T)-  IK?.  Kj)  , 


-  K^(Kj)  . 


^  ^  «f> 


(63) 


which  are  represented  in  Figs*  2  through  4  and  are  valid  in  the  interval 
0  <  <  9A0#  Incidentally,  the  solution  corresponding  to  ■  1/2  is  a 

wedge  • 

Bodies  of  class  II >  These  bodies  consist  of  a  reguleu:  shape  followed 
by  a  constant  thickness  contour  and  are  obtained  for  0  ^  5^  <  1  and  K2  * 

The  shape  of  the  optimum  body,  the  thickness  parameter,  the  drag  coefficient  * 
and  the  friction  parameter  are  written  as 


0  s  ?  i  ?^  . 


71-1 


(64) 


Elimination  of  the  abscissa  of  the  transition  point  from  these  equations 
leads  once  more  to  functional  relationships  of  the  form  (63)  which  are  plotted 
In  Figs*  2  through  4  and  are  valid  for  9/10  <  As  the  friction  pa¬ 

rameter  increases^  the  abscissa  of  the  tremsitlon  point  moves  forwardf  s  clr- 
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oufflstanot  which  la  common  to  all  of  the  boundary  value  problama  analysed 
here# 


6#5«  given  Thickneaa  and  ^k)ment  of  Inertia  of  the  Contour 

If  the  thickneaa  and  the  moment  of  inertia  of  the  contour  are  preacribed 
while  the  encloaed  area  la  freei  the  tranaveraality  condition  leads  to  K2  ■  0« 
Aa  in  the  previous  case^  a  one->paraineter  family  of  solutions  exists#  These 
solutions  are  of  thither  class  1  or  class  II  depending  on  whether  the  friction 
coefficient  Is  smaller  or  larger  than  a  certain  critical  value#  Once  more, 
the  representation  of  the  results  is  facilitated  by  Introducing  a  thickness 
parameter  and  a  friction  parameter.  These  parameters  are  defined  as 

Kt  ■  ^  T  , 

and,  because  of  £qs#  (35)  and  (38-2),  can  be  rewritten  as 

Bodies  of  Class  I#  These  bodies  consist  of  a  regular  shape  only  and 
are  obtained  for  ^  and  -  «  <  1#  Using  £q#  (32-1) «  one  can  express 

the  geometry  of  the  optimum  shape  In  the  form 


(67) 


where 


g(iii  K3) 


g(71,  Kj)  -  * 


(68) 
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^  r(<p,  k)  -  ^  B(<p.  k) 


^5“ T 1  ± y  1  - KjTi^ 


and  where  the  upper  algn  is  valid  for  ^  0  and  the  lower  ai^i  for  2:  0* 
Ihe  aymbols  F  and  £  denote  the  incomplete  elliptic  integrals  of  the  first 
and  second  kind  whose  argtuaents  tp  and  parameter  k  are  defined  as 


Kj) 


«  arc  cos 


A 1 T  y  1  - 

1  ±  y  1  -  KjTi^ 


(69) 


Because  of  Eqs*  (46)  and  (65)t  the  thickness  parameteri  the  drag  coefficient^ 
and  the  friction  parameter  can  be  expressed  in  the  fora 


>  7K, 


(1  -  Kj) 

^  3g(i,  kj)  ■ 


2/y 


7  ■  §  ^5^  IS>  ♦ 


(70) 


Elimination  of  the  paurameter  from  Eqs*  (6?)  and  (70)  yields  functional 
relationships  of  the  form  (63)  which  are  plotted  in  Figs«  5  through  7  and  are 
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valid  in  the  Interval  0  <  <  0»33*  Incldentallj,  the  solution  correspond* 

ing  to  «  1/3  is  a  wedge* 

Bodies  of  olaee  II*  These  bodies  consist  of  a  regular  shape  followed 
by  a  constant  thickness  contour  and  are  obtained  for  0  <  <  1  and  ■  1* 

After  the  shape  of  the  optinun  body,  the  thickness  parameter^  the  drag  coef¬ 
ficient!  and  the  friction  parameter  aure  written  as 

0  <  ?  <  5c  .  ^ 

5c  <  ?  *  1  f  T! 


“5 


elimination  of  the  abscissa  of  the  transition  point  from  these  equations 
leads  once  more  to  functional  relationships  of  the  form  (63)  which  are 
plotted  in  Figs*  5  through  7  ud  are  valid  in  the  interval  0*55  ^  ^ 

6*6*  Given  Enclosed  Area  and  Moment  of  Inertia  of  the  Contoiir 

If  the  enclosed  area  and  the  moment  of  inertia  of  the  contour  are  given 
while  the  thickness  is  fresi  the  transversal! ty  condition  leads  to 


(71) 


Ml 
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«  1  -  K^«  Again,  tha  extremal  arc  is  represented  by  a  one-parameter 
family  of  aolutione  of  either  class  I  or  class  II  depending  on  whether  the 
friction  coefficient  Is  smaller  or  larger  than  a  certain  critical  value# 

Once  more,  the  representation  of  the  results  is  facilitated  if  a  thickiiess 
parameter  and  a  friction  parameter  are  introduced#  These  parameters  are  de¬ 
fined  by 

K  -  T  ,  K.  -  (72) 


and,  because  of  Eqs#  (38) ,  can  be  rewritten  as 


WV 


Bodies  of  class  I#  These  bodies  consist  of  a  regular  shape  only  and 
are  obtained  for  5^  “  ^  »>Kji-l#  If  Eq#  (32-1)  is  employed,  the 
expression  for  the  optimum  shape  can  be  written  as 


h(0,  K3)  ^  h(T|,  K3) 

?  -  RoTk^TTWT^) 


(74) 


In  the  numerator  of  the  above  equation,  the  upper  sign  is  V8j.id  for 
T)  *  (Kj  -  1)/2Kj,  and  the  lower  sign,,  for  T1  i  (K^  -  l)/2Kj  if  •  2  1^  k  1{ 
for  1  2  2  -  1,  the  upper  sign  is  valid  for  all  values  of  T|,  With  regard 

to  the  denominator,  the  lower  sign  is  valid  when  «  2  K,  2  1,  and  the  upper 


30 


sign,  when  1  %  -  1.  The  defixiition  of  the  fwetlon  h(11y  K^)  which 

appears  in  Eq.  (7^)  is  given  by 


Mu,  Kj) 


^  r(q),  k)  -  ^  B(q),  k) 

1  -  a 


(75) 


where  the  upper  sign  is  valid  for  «  %  >  !«  the  lower  sign  is  valid 

for  1  >  K,  >  -  1,  and  the  quantities  or,  (p,  and  k  are  defined  as 


or(71,  Kj) 


~4k^(1  -t-  K3T))  (1  -  Tp 

.  (1  +  Kj)^ 


(76) 


MTl.  Kj) 


«  arc  COB 


V?  +  1  -  Of 


f 


Finallj^  because  of  Eqs*  (46)  and  (73)  f  the  thickness  parameter^  the  drag 
coefficient^  and  the  friction  parameter  can  be  expressed  in  the  form 


(77) 


where 
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•<ICj)  -  3  [wo.  tj)  *  Wl.  *5>] 

'<«3>  ■  ^  l- « -  ■s’C""’*  IS’  ’  vij  * 

“■b*  ■  7^  ['  - 1<‘  -  *5>  ’] 


(7«) 


and  whara  the  lower  signs  are  valid  for  «  >  >  1  and  tha  upper  signs^ 

for  1  >  Im  Eliisination  of  the  parameter  from  Eqs*  (7^)  and  (77) 

yields  functional  relationships  of  the  form  (63)  which  are  plotted  in  figs*  8 
through  10  and  eure  valid  for  0  < 

Bodies  of  class  II#  These  bodies  consist  of  a  regular  shape  followed 
by  a  constant  thickness  contour  and  are  obtained  for  0  <  <  It  K2  * 

and  »  «  !•  After  the  shape  of  the  optimum  bodyt  the  thickness  parameter , 
the  drag  coefficient t  and  the  friction  parameter  are  written  as 

,  mJIiVL 

(79) 


147  ~ 

^  (14  - 
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elimination  of  the  absciesa  of  the  transition  point  from  these  equations 
leads  once  more  to  functional  relationships  of  the  form  (63)  which  are 
plotted  in  Figs*  8  through  10  and  are  valid  for  ^9/16  ic  •• 
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7#  DISCUSSION  AND  CONCLUSIONS 

From  the  previous  analyslsi  it  appears  that,  despite  the  gweralitj 
of  the  present  problem,  the  method  of  solution  is  relatively  simple  and 
has  the  merit  of  leading  to  analytical  solutions  in  each  of  the  six  par¬ 
ticular  cases  considered  here*  The  main  comments  to  these  solutions  are  as 
follows: 

(a)  For  the  general  problem  in  which  the  length  is  free  and  arhitrarj 
conditions  are  assigned  to  the  thickness,  the  enclosed  area,  and  the  moment 
of  inertia  of  the  contour,  the  totality  of  extremal  arcs  is  represented  by 

a  two-parameter  family  of  solutions  if  dimensionless  coordinates  are  employed, 
that  is,  if  the  abscissa  and  the  ordinate  are  normalized  with  respect  to 
the  length  and  the  semi  thickness*  Each  member  of  the  family  is  characterized 
by  a  sharp  leading  edge*  Furthermore,  each  extremal  arc  may  involve  at  most  one 
corner  point  and,  hence,  two  subarcs*  Of  these  subarcs,  one  is  characterized 
by  a  positive  pressure  coefficient  and  is  called  the  regular  shape;  the  other 
is  characterized  by  a  zero  pressure  coefficient  and  is  called  the  zero-slope 
shape*  Consequently,  two  classes  of  bodies  can  be  identified:  (I)  bodies 
composed  of  a  regular  shape  only  and  (II)  bodies  composed  of  a  regular  shape 
followed  by  a  constant  thickness  contour* 

(b)  If  only  one  geometric  quantity  is  assigned  (the  thickness,  the 
enclosed  area,  or  the  moment  of  inertia  of  the  contour),  a  zero-parameter 
family  of  solutions  exists  (that  is,  a  single  curve)*  In  all  cases,  the  so¬ 
lution  is  of  class  I,  that  is,  consists  of  a  regular  shape  only*  In  particu¬ 
lar,  if  the  thickness  is  given,  the  solution  is  a  wedge,  and  its  length  is 
such  that  the  friction  drag  is  2/3  of  the  total  drag*  If  the  enclosed  area 
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ie  given  I  the  coaplemente  of  the  ordinate  and  the  absoisea  obey  a  3/2  power 
laW|  and  the  optlaum  thickness  ratio  is  such  that  the  friction  drag  is  5/^ 
of  the  total  drag.  Finally |  if  the  moment  of  Inertia  of  the  contour  ie 
given t  the  solution  is  represented  by  a  combination  of  elliptic  integrals 
of  the  first  and  the  second  kind|  and  the  optimum  thickness  ratio  is  such 
that  the  friction  drag  is  7/9  of  the  total  dreig* 

(c)  If  two  geometric  quantities  eire  prescribed  (the  thickness  and  the 
enclosed  aresi  the  thickness  and  the  moment  of  inertia  of  the  contour ^  or 
the  enclosed  area  and  the  moment  of  inertia  of  the  contour),  a  one-paraneter 
family  of  solutions  exists.  This  parameter,  called  the  friction  parameter, 
ie  proportional  to  the  cubic  root  of  the  friction  coefficient  and  is  indicative 
of  the  relative  Importance  of  the  friction  drag  with  respect  to  the  pressure 
drago  Depending  on  the  value  of  the  friction  parameter,  two  distinct  behaviors 
are  possible.  If  the  friction  parameter  is  subcritical  (smaller  than  a  certain 
critical  value),  the  solution  is  of  class  I  and,  therefore,  involves  a  regu¬ 
lar  shape  only.  If  the  friction  parameter  is  supercritical  (larger  than  a 
certain  critical  value),  the  solution  is  of  class  II  and,  therefore,  in¬ 
volves  a  regular  shape  followed  by  a  constant  thickness  contour;  in  all 
cases,  the  trr^)  ^ition  point  from  the  regular  shape  to  the  constant  thickness 
contour  shifts  forward  as  the  friction  parameter  increases. 

In  closing,  it  should  be  noted  that,  if  the  limiting  process  0  is 
carried  out,  the  present  solutions  reduce  to  the  inviscid  flow  solutions 
already  calculated  in  Ref.  1*  It  should  also  be  noted  that  some  of  the  opti¬ 
mum  shapes  obtained  with  this  analysis  aure  concave;  consequently,  these 
bodies  should  be  restudied  using  the  Newton-Busemann  pressure  coefficient 
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law;  this,  howevery  requires  a  more  thorough  understanding  of  the  friction 
drag  associated  with  the  possible  presence  of  free  layers*  Finallyi  when 
the  square  of  the  thickness  ratio  becomes  nonnegllgible  with  respect  to 
one,  the  slender  body  approximation  is  violated;  consequently,  the  problem 
should  be  reinvestigated  using  the  exact  Newtonian  expression  for  the  pres-- 
sure  coefficient,  that  is,  the  sine  square  law* 
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Fig.  1.  Optimum  shapes  for  given  values  of  the  thickness,  the  enclosed 
area,  or  the  moment  of  inertia  of  the  contour. 


Fig.  2.  Optimum  shapes  for  given  thickness  and  enclosed  area. 


Fig.  3.  Thickness  ratio  for  given  thickness  and  enclosed  area. 
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Fig.  6.  Thickness  ratio  for  given  thickness  and  moment  of  inertia  of  the 
contour. 
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Fig.  7.  Drag  coefficient  for  given  thickness  and  moment  of  inertia  of  the 
contour. 


Fig.  8.  Optimum  shapes  for  given  enclosed  area  and  moment  of  inertia  of 
the  contour. 


Fig.  9 .  Thickness  ratio  for  given  enclosed  area  and  moment  of  inertia  of 
the  contour. 
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Fig.  10.  Drag  coefficient  for  given  enclosed  area  and  moment  of  inertia  of 
the  contour. 


